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ABSTRACT 

The origin and initial mass function of young stellar clusters near the Galactic 
center are still poorly understood. Two of the more prominent ones, the Arches 
and Quintuplet clusters, may have formed from a shock-induced burst of star for- 
mation, given their similar age and proximity to each other. Their unusual mass 
distribution, however, may be evidence of a contributing role played by other fac- 
tors, such as stellar capture from regions outside the clusters themselves. Diffuse 
X-ray emission from these sources provides us with a valuable, albeit indirect, 
measure of the stellar mass-loss rate from their constituents. Using recent data 
acquired with Chandra, we can study the nature and properties of the outflow 
to not only probe the pertinent physical conditions, such as high metallicity, the 
magnetic field, and so forth, but also to better constrain the stellar distribution 
within the clusters, in order to identify their formative history. In this paper, 
we present a set of three-dimensional smoothed particle hydrodynamics simula- 
tions of the wind-wind interactions in both the Arches and Quintuplet clusters. 
We are guided primarily by the currently known properties of the constituent 
stars, though we vary the mass-loss rates in order to ascertain the dependence 
of the measured X-ray flux on the assumed stellar characteristics. Our results 
are compared with the latest observations of the Arches cluster. Our analysis of 
the Quintuplet cluster may be used as a basis for comparison with future X-ray 
observations of this source. 

Subject headings: galaxies: clusters: individual (Arches, Quintuplet) — Galaxy: 
center — radiation mechanisms: thermal — shock waves — stars: winds — X-rays: 
diffuse 
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1. Introduction 

Understanding the environment's role in star formation and, in turn, the feedback ex- 
erted by star formation on the Galactic environment, is a problem of significance to several 
fields in astronomy, from the creation of compact objects to the formation of galaxies. The 
Galactic center, with its relatively high magnetic field strength, clouds with high particle 
density, and large velocity dispersions, provides an ideal environment to study star formation 
under extreme conditions (Morris 1993; Melia & Falcke 2001). This type of investigation 
can benefit from the existence of several stellar clusters in this region, including the Arches 
and Quintuplet clusters. Learning more about their stellar constituents, and possibly their 
formative history, may even give us a glimpse into the emergence of objects that will ulti- 
mately populate the central parsec of the Galaxy (e.g., Gerhard 2001; McMillan & Portegies 
Zwart 2003). 

The Arches and Quintuplet clusters have been studied over a range of wavelengths, 
from radio to X-ray. X-rays provide a unique window for investigating both the formation 
of binaries (point sources from binary interactions) and the wind interactions within the 
entire cluster (diffuse emission). As Rockefeller et al. (2004) have shown, the diffuse X-ray 
emission is a sensitive measure of the mass-loss rate of stars in mutually interactive situations. 
Stellar mass-loss remains one of the largest uncertainties in stellar evolution (e.g., Wellstein 
& Langer 1999). X-ray observations of these clusters represent a unique probative tool for 
studying the winds produced by high-metallicity systems. 

In this paper, we model the propagation and interaction of winds from stars in both 
the Arches and Quintuplet clusters, calculating the X-ray fluxes arising from the consequent 
shocked gas. Yusef-Zadeh et al. (2002) serendipitously discovered the Arches cluster with 
Chandra and identified five components of X-ray emission, which they labeled A1-A5, though 
only A1-A3 seem to be directly associated with the cluster (see Figure 2, Yusef-Zadeh et al. 
2002). Al is apparently associated with the core of the cluster, while A2 is located ~ 10" 
northwest of Al. Al and A2 are partially resolved, while A3 is a diffuse component that 
extends beyond the boundary of the cluster; Yusef-Zadeh et al. (2002) speculated that some 
or all of A1-A3 may be produced by interactions of winds from stars in the system. Analysis 
of additional Chandra observations that covered the Arches cluster and first results from the 
Quintuplet cluster are presented by Wang (2003) and Law & Yusef-Zadeh (2004); these new 
observations resolve Al into two distinct components, labeled AlN and AlS, and indicate 
that AlN, AlS, and A2 are all point-like X-ray sources. 

Several efforts have already been made to study the X-ray emission from clusters. Oz- 
ernoy et al. (1997) and Canto et al. (2000) performed analytic calculations to estimate the 
diffuse emission from these systems. The interaction of winds in the Arches cluster has 
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been simulated by Raga et al. (2001) using the "yguazii-a" adaptive grid code described 
in Raga et al. (2000). However, all previous work focused exclusively on the Arches clus- 
ter, and even the detailed simulations assumed identical large values for the mass-loss rates 
(M = 10~^Mq yr~-^) and wind velocities {v^ind = Ij 500 km s^^) of the constituent stars. In 
this paper, we present the results of simulations of both the Arches and Quintuplet clusters, 
using detailed radio flux measurements (where available) and spectral classifications to pin 
down the expected mass-loss rates of stars in each system. We then compare our results to 
the most recent X-ray observations of these two clusters, including new data presented in 
this paper. 

A summary of relevant properties of each cluster is presented below. We describe our 
numerical technique, including the characteristics of the wind sources and the gravitational 
potential of the clusters, in § 2 of the paper. The new observations are described in § 3. A 
comparison of the theoretical results with the data is made in § 4, and the relevance to the 
Galactic center conditions is discussed in § 5. 

1.1. The Arches Cluster 

The Arches stellar cluster is an exceptionally dense aggregate of stars located at I — 
0.12°, b — 0.02°, about 11' in projection from the Galactic center (see e.g. Nagata et al. 1995; 
Cotera et al. 1996; Figer et al. 2002). The cluster is apparently a site of recent massive star 
formation; it contains numerous young emission-line stars which show evidence of strong 
stellar winds. Using near-IR color-magnitudes and K band counts, Serabyn et al. (1998) 
estimated that at least 100 cluster members are O stars with masses greater than 20 Mq 
and calculated a total cluster mass of ~ (1.5-6) x 10^ Mq. Figer et al. (1999b) used HST 
NICMOS observations to determine the slope of the initial mass function (IMF) of the Arches 
cluster and calculated a cluster mass of ~ 10^ M©, with possibly 160 O stars and an average 
mass density of ~ 3 x 10^ Mq pc~^. 

The 14 brightest stars of this cluster have been identified with JHK photometry and 
Bra and Br7 hydrogen recombination lines, showing that they have the characteristic colors 
and emission lines of Of-type or Wolf-Rayet (WR) and He I emission-hne stars. Nagata et al. 
(1995) inferred from the strength of the Bra and Br7 hne fluxes that these stars are losing 
mass at a prodigious rate, M ~ 2 x 10~^ Mq yr~^, in winds moving at ~ 10^ km s~^. Cotera 
et al. (1996) confirmed the presence of young, massive stars using K-hand spectroscopy; they 
identified 12 stars in the chistcr with spectra consistent with late-type WN/Of objects, with 
mass-loss rates M ~ (1-20) x 10~^ Mq yr~^ and wind velocities v^o ~ 800-1, 200 km s~^. 
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FoUow-up Very Large Array (VLA) observations at centimeter wavelengths of the bright- 
est stars in the cluster have solidified the detection of powerful ionized stellar winds. Using 
the observed 8.5 GHz fiux densities from 8 sources in the Arches cluster and the relationship 
between flux density and mass- loss rate derived by Panagia & Felli (1975), 

M = (5.9 X 10-^ Mo yr-i) ( ( ( -^—\ (1) 

^ ^VlmJyy/ V500kms-V Vskpcy/ ' ^' 

where 5*8.5 is the 8.5 GHz fiux density, Voo is the wind terminal velocity, and d is the distance 
to the source (~ 8 kpc, for stars in the Arches cluster), Lang et al. (2001) calculated mass 
loss rates M — (3-17) x 10~^ Mq yr~^, assuming a wind electron temperature T ~ 10^ K, 
Z — 1, and a mean molecular weight /j, — 2. The Wolf-Rayet phase is short-lived, but while 
in this mode, stars dominate the mass ejection within the cluster. 

In an attempt to represent both the identified stellar wind sources and the population 
of stars likely to be producing significant but currently undetected winds, we include 42 
wind sources (listed in Table 1 and shown in Figure 1) in our simulations of the Arches 
cluster. The stars labeled AR1-AR9 correspond to the 9 radio sources identified by Lang 

(2002). The first 29 stars (labeled 1-29) in Table 1 have mass estimates greater than 60 Mq 
(Figcr ct al. 2002), and are likely the most powerful sources of wind in the cluster. The 
remaining 13 stars used in the simulations have masses less than 60 Mq but greater than 
25 Mq, and are located on the north side of the cluster; they arc included to better represent 
the spatial pattern of X-ray emission around the core of the cluster. Based on the broadening 
of the Br7 fine observed by Cotera et al. (1996), stars in the Arches cluster simulations are 
assigned wind velocities of 1, 000 km s~^. The stars labeled AR1-AR9, which have observed 
8.5 GHz fiux densities, are assigned mass- loss rates according to Equation 1. Stars with no 
associated 8.5 GHz detection but with masses larger than 60 Mq are assigned a mass-loss 
rate of 3 x 10~^ Mq yr~^, which is equal to the lowest mass-loss rate inferred from the 
weakest observed 8.5 GHz signal from the Arches cluster (Lang 2002). Stars with masses 
less than 60 Mq are assigned a mass-loss rate of 3 x 10^^ Mq yr^^; their winds will have 
little effect on the overall luminosity but may alter the shape of the X-ray-emitting region. 



1.2. The Quintuplet Cluster 

Slightly further north of Sgr A*, the Quintuplet cluster is located at / = 0.16°, b = 0.06°. 
The cluster has a total estimated mass of --^ 10^ Mq and a mass density of ~ 10^ Mq pc^^ 
(Figer et al. 1999a). Like the objects in the Arches cluster, the known massive stars in the 
Quintuplet cluster have near-IR emission-line spectra indicating that they too have evolved 
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away from the zero-age main sequence and now produce high-velocity stellar winds with 
terminal speeds of 500-1, 000 km s~^. 

Figer et al. (1999a) obtained K-band spectra of 37 massive stars in the Quintuplet 
cluster and found that 33 could be classified as WR stars, OB supergiants, or luminous blue 
variables (LBVs), implying a range of wind mass-loss rates M ~ (0.1-6.6) x 10~^ Mq yr~^. 
VLA continuum images at 6 cm and 3.6 cm of the Sickle and Pistol H ll regions reveal eight 
point sources located in the vicinity of the Quintuplet cluster, including the radio source at 
the position of the Pistol nebula (Lang et al. 1999). These are labeled QRl through 7, and 
the Pistol star, in Figure 2 below. 

The near-IR counterparts of QR4 and QR5 are hot, massive stars with high mass-loss 
rates, one an OB I supergiant and the other a WN9/0fpe (Figer et al. 1999a). The sources 
QRl, QR2, and QR3 also have spectral indices consistent with stellar wind sources, but they 
have no obvious NICMOS stellar counterparts. Lang et al. (1999) speculated that this may 
be due to variable extinction associated with a dense molecular cloud located in front of 
the cluster. Given the uncertainty, and possible variation, in stellar identification for these 
5 objects, we therefore adopt a value of ~ 500 km (typical in OB supergiants) for the 
speed of their wind. The Pistol star, on the other hand, is a prominent source in the near-IR 
NICMOS image, and is evidently a luminous blue variable (Figer et al. 1998) possessing a 
powerful stellar wind, though with a terminal speed of only v^o ~ 100 km s~^. 

In our simulations of the Quintuplet cluster, we include 31 massive stars with spectral 
classifications identified by Figer et al. (1999a), using estimates for the wind parameters of 
each star according to its spectral type. Wind velocities are determined according to the 
broad spectral type of each star: we assume a velocity of 1, 000 km s^^ for winds from WR 
stars, 500 km s~^ winds for OB supergiants, and 100 km s~^ for the Pistol star, a LBV. 

For those stars that are radio sources (QR1-QR3, QR6, and QR7 from Lang et al. 
1999), mass- loss rates are determined according to Equation 1. For the rest of the stars, 
we estimate the mass- loss based on the spectral classifications by Figer et al. (1999a). For 
OB stars, the mass-loss rate was assumed to be 10^^ Mq yr^^ for stars with classification 
higher than BO. For smaller stars, the mass-loss rate was assumed to be 10~® Mq yr~^. 
For Wolf-Rayet stars, we use a luminosity-mediated mass-loss relation (Wellstein & Langer 
1999): 

where L is the stellar luminosity from Figer et al. (1999a), Xg is the hydrogen mass fraction, 
and A; is a constant which we calibrated using our radio-determined mass-loss rates. The 
locations and wind parameters of the stellar wind sources are summarized in Table 2, and 
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their relative positions are shown in Figure 2. 



2. The Physical Setup 

Our calculations use the three-dimensional smoothed particle hydrodynamics (SPH) 
code discussed in Fryer & Warren (2002) and Warren et al. (2004), modified as described 
in Rockefeller et al. (2004) to include stellar wind sources. The gridless Lagrangean nature 
of SPH allows us to concentrate spatial resolution near shocks and model gas dynamics and 
wind-wind interactions on length scales that vary by several orders of magnitude within a 
single calculation. 

We assume that the gas behaves as an ideal gas, according to a gamma- law (7 = 5/3) 
equation of state. The effect of self-gravity on the dynamics of the gas should be negligible 
compared to the effect of the central cluster potential; we calculate gravitational effects by 
approximating the potential of each cluster with a Plummer model. 

The computational domain for the simulation of each cluster is a cube approximately 
6 pc on a side, centered on the middle of the cluster. To simulate "flow-out" conditions, 
particles passing through the outer boundary are removed from the simulation. The initial 
conditions assume that the space around and within each cluster is empty; massive stars 
in each cluster then inject matter into the volume of solution via winds as the calculation 
progresses. The number of particles in each simulation initially grows rapidly, but reaches 
a steady number (~ 6.6 million particles for simulations of the Arches cluster and ~ 3.3 
million particles for the Quintuplet simulations, since there are fewer identified wind sources 
in the latter) when the addition of particles from wind sources is compensated by the loss of 
particles through the outer boundary of the computational domain. 



2.1. Cluster Potential 

We model the gravitational potential $ of a cluster using a Plummer model (Plummer 
1911), 

-CM 

where M is the total mass of the cluster. The radial density profile assumed for the cluster 
is therefore 

/ 2 \ -5/2 

p(r)=Po(l + ^) , (4) 
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where the central density po is 

We note that the spatial distributions of massive stars in the clusters are not entirely con- 
sistent with the distribution implied by the Plummer potential; for example, the average 
projected distance of both the set of massive stars and the set of all observed stars in the 
Arches cluster is roughly twice the average distance predicted by a Plummer model based 
on the estimated total mass and central density of the cluster. We include the Plummer po- 
tential to approximate the combined gravitational influence of the entire cluster, including 
the estimated mass of stars too small or dim to be observed. 

Figer et al. (1999a) provide estimates of the total mass of each cluster by measuring the 
mass of observed stars, assuming a Salpeter IMF slope, and extrapolating down to 1 Mq — 
observed stars account for at most 25% of the mass of the Arches cluster and 16% of the mass 
of the Quintuplet cluster. They also estimate the density of each cluster by determining the 
volume of the cluster from the average projected distance of stars from the cluster center, 
and dividing the total mass by this volume; because the values are calculated using the total 
cluster mass but only the average projected cluster radius, the density estimates are probably 
closer to the central densities than the average densities. We assume that the values reported 
are the central densities of the clusters. The total cluster mass and central density and the 
calculated value of b for each cluster are presented in Table 3. 



2.2. Wind Sources 

We implement wind sources as literal sources of SPH particles, using the scheme de- 
scribed in Rockefeller et al. (2004). The mass loss rates and wind velocities inferred from 
observations are reported in Tables 1 and 2. We position each source at its observed x and 
y location and choose the z coordinate randomly, subject to the constraint that the wind 
sources are uniformly distributed over a range in z equal to the observed range in x and y. 
The choice of z positions has a much smaller effect on the X-ray luminosity than the choice 
of mass-loss rate (discussed below); Rockefeller et al. (2004) performed simulations of wind 
sources in the central few parsecs of the Galaxy with two different sets of z positions, and 
the average 2 — 10 keV X-ray luminosity from the central 10" of the simulations differed by 
only 16% (7.50 x 10^^ erg s~^ arcsec"^ from a simulation with a dense arrangement of wind 
sources in the center of the volume of solution versus 6.45 x 10^^ erg s'^ arcsec"^ from a 
simulation with wind sources uniformly distributed in z). In addition, Raga et al. (2001) 
conducted three simulations of the Arches cluster in which z positions of sources were chosen 
by sampling from a distribution function f{R) oc R~^; they found a difference in 0.5 — 8 keV 
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X-ray luminosity of only 3% between the most and least luminous simulations. 

Figure 1 shows the positions in the sky plane of the wind sources in the Arches cluster, 
while Figure 2 shows the wind sources in the Quintuplet cluster; the size of the circle marking 
each source corresponds to the relative mass loss rate (on a linear scale) for that star. The 
initial temperature of the winds is not well known; for simplicity, we assume that all of the 
winds have a temperature of 10^ K. Our results arc fairly insensitive to this value, however, 
since the temperature of the shocked gas is determined primarily by the kinetic energy flux 
carried into the collision by the winds. The sources are assumed to be stationary over the 
duration of the simulation. 

We conduct two simulations of the wind-wind interactions in each cluster; the simula- 
tions differ only in the choice of mass-loss rate for each star, with the wind speed assumed 
constant. The "standard" simulations use mass- loss rates inferred from observations as de- 
scribed in §§ 1.1 and 1.2, while the "high-M" simulations use mass-loss rates increased from 
the standard value by a factor 2. 

3. Observations 

An on-axis Chandra observation of the Arches cluster was taken on 2004 June 8 for 
98.6 ksec. The ACIS-I detector was placed at the focal plane. The data were analyzed 
with the latest CIAO (version 3.1). While results of this observation will be presented in an 
upcoming paper (Wang 2004), we concentrate here on the comparison of the data with the 
simulations. Figure 3 shows an overlay of an X-ray intensity contour map on a HST NICMOS 
image of the Arches cluster (Figer et al. 1999b). The actual spatial resolution (with a FWHM 
< 1") is better than what appears on this smoothed image. The position coincidence of 
the point-like X-ray sources and bright near-IR objects is apparent. To compare with the 
simulated cluster wind properties, we remove a region of twice the 90% energy-encircled 
radius around each of the sources. 



4. Results 

All four simulations were run significantly past the point in time when the stellar winds 
fill the volume of solution and gas shocked in wind- wind collisions fills the core of each cluster. 
The Arches cluster simulations were run past t = 10, 000 yr; the Quintuplet simulations were 
run for over 14, 000 yr. The winds in each simulation reach the edge of the volume of solution 
after ~ 3, 000 yr, but the most relevant timescale for determining when the simulations reach 
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a steady level of X-ray luminosity is the time required to fill the core of each cluster with 
shocked gas. The Arches cluster core is roughly five times smaller in radius than the core of 
the Quintuplet cluster, so the X-ray luminosity from the Arches simulations reaches a steady 
value relatively quickly compared to the Quintuplet simulations. 

4.1. Total Flux and Time Variation 

In order to calculate the observed continuum spectrum, we assume that the observer is 
positioned along the positive 2;-axis at infinity and we sum the emission from all of the gas 
injected by the stellar wind sources into the volume of each calculation. For the conditions 
encountered in the two clusters, scattering is negligible and the optical depth is always less 
than unity. For these temperatures and densities, the dominant components of the continuum 
emissivity are electron-ion (egj) and electron-electron (eee) bremsstrahlung. We assume that 
the gas is in ionization equilibrium, although spectra obtained from the Arches cluster with 
Chandra indicate the presence of line emission at 6.4 keV; however, no other significant lines 
are observed, so our use of a bremsstrahlung model without line emission is a reasonable 
approximation. Future, more refined versions of the calculations reported here will need to 
include the effects of partial ionization in the shocked gas, a condition suggested by the iron 
line emission. We note, however, that the overall energetics and diffuse X-ray luminosity 
calculated by assuming full ionization err only marginally since the fraction of charges free 
to radiate via bremsstrahlung should be very close to one. 

The X-ray luminosities calculated from our simulations support the recent result of 
Law & Yusef-Zadeh (2004) that the majority of the X-ray emission from these clusters (e.g., 
~ 60% for Arches) is probably due to point sources. The diffuse 0.2-10 keV X-ray luminosity 
calculated from our "standard" Arches cluster simulation — 5.4 x 10^^ erg s~^ — falls below 
the 0.2-10 keV luminosity of 4.1 x 10^^ erg s~^ from the Al and A2 components reported 
in Yusef-Zadeh et al. (2002); the simulation using elevated estimates for the mass- loss rates 
produces 2.2 x 10^^ erg between 0.2 and 10 keV, or 53% of the emission observed by 
Chandra. On the other hand. Law Sz Yusef-Zadeh (2004) now identify Al and A2 as point- 
like components; after subtracting the contributions of Al and A2, Yusef-Zadeh et al. (2002) 
find that the A3 component has a 0.5-10 keV luminosity of ~ 1.6 x 10'^^ erg s^^. Our 
"standard" simulation produces 2.7 x 10'^"^ erg s^^ between 0.5 and 10 kcV; lowering the 
mass-loss estimates of all wind sources by ~ 30% or slightly decreasing the assumed wind 
velocities would produce even better agreement. 

Our "standard" and "high-M" simulations of the Quintuplet cluster produce 1.5 x 
-|^g33 gj,g g-1 5 g ^ ^q33 gj,g respectively, between 0.5 and 8 keV. Law & Yusef-Zadeh 
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(2004) estimate a 0.5-8 keV luminosity of ~ 1 x 10^^ erg s~^ for the diffuse emission from 
the Quintuplet cluster, but they point out that other regions of diffuse emission to the north 
and south of the cluster introduce a complicated pattern of background emission and limit 
the precision of this estimate. 

The long-term variation of the X-ray luminosity as a function of time demonstrates some 
of the key differences between the two clusters. Figures 4 and 5 show the time variation of 
the 0.5-8 keV X-ray luminosity from all four simulated clusters. Each large graph shows 
the variation over the course of the entire simulation (> 10, 000 yr for the simulations of 
the Arches cluster, and > 14,000 yr for the Quintuplet simulations), while each inset graph 
shows variation over 1,000 timesteps, or ~ 150 yr. The fact that the Quintuplet cluster 
is nearly five times larger in radius than the Arches cluster means that much more time is 
required for shocked gas to fill the central region. The Arches cluster reaches a steady X-ray 
luminosity after less than 2, 000 yr, while the luminosity of the Quintuplet cluster does not 
clearly level off until more than 10, 000 yr have passed. 

The two clusters also exhibit differences in the size of short-term fluctuations in X-ray 
luminosity. Both simulations of the Arches cluster show short-term variations in luminosity 
of ~ 1% over timescales of ~ 50 yr, while the X-ray luminosity in the "standard" Quintuplet 
simulation fluctuates by ~ 4%, and short-term variations in the "high-M" simulation of the 
Quintuplet cluster are as large as ~ 7%. 

4.2. Spatial Variation of the X-ray Flux 

Figure 3 shows that many of the bright X-ray peaks in the Arches cluster correspond to 
actual stars, presumably binaries whose binary wind interactions produce the strong localized 
X-ray emission. These sources must be subtracted to study the diffuse X-ray emission. 

The simulated 0.5-8 keV X-ray contours from the region near the core of the Arches 
cluster, shown in Figure 6, are generally comparable to the contours generated from Chandra 
observations (Figure 3; see also Figure 2b, Yusef-Zadeh et al. 2002). The strongest emission 
in the simulations and in Chandra images comes from the core of the cluster, and both 
sets of contours form elliptical patterns aligned primarily north-south. X-ray contours from 
the simulation of the Quintuplet cluster are shown in Figure 7; the Quintuplet cluster is 
significantly less dense than the Arches cluster, so the X-ray emission is correspondingly less 
strongly peaked toward the center of the cluster. 

The plots in Figure 8 show the total 0.5-8 keV luminosity from concentric columns 
aligned along the line of sight and extending outward in radius from the center of the Arches 
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(left) and Quintuplet (right) clusters. The most luminous gas in the Arches cluster is confined 
to within 20" of the center of the cluster; in contrast, we include stars beyond 50" in the 
Quintuplet cluster, and its luminosity continues to increase even beyond a radius of 50". 

Similarly, the plots in Figure 9 show the 2-8 keV X-ray flux per square arcmin as a 
function of distance from the center of each cluster from all four simulations. The individual 

crosses and error bars in the graph from the Arches cluster represent flux measurements from 
Chandra observations of the cluster, after point sources have been removed and an estimate 
of the background X-ray flux has been subtracted. Here we assume that all of the emission 
centered around NICMOS stellar sources is due to binary wind interactions. The estimated 
background — 0.064 counts s^^ arcmin"^ — is the average number of counts obtained in an 
annulus between radii of 0.5' and 0.9'. The simulations of the Arches cluster apparently 
produce more X-rays near the center of the cluster but decrease in intensity more rapidly 
toward larger radii. The relatively flat surface brightness proflle evident in the Chandra data 
outside a radius of 0.3' may arise in part from the presence of additional background X-ray 
emission near the Galactic center. It may also indicate confinement of the X-ray-emitting 
gas in the Arches cluster by ram pressure exerted by a molecular cloud or other external 
medium surrounding the cluster; our calculations include no such medium, so the simulated 
gas escapes and cools more rapidly as it leaves the core of the cluster. 

One way to constrain the amount of confining material surrounding the cluster is to 
measure expansion velocities of the gas out of the cluster. In our simulations, which did 
not include confining gas, the material accelerates until it reaches an asymptotic velocity 
limit roughly equal to the mean wind velocity (Figure 10). In the Arches cluster, we reach 
that limit. In the (larger) Quintuplet cluster, that limit apparently occurs at a distance 
from the cluster center that is larger than the size of the simulation. However, if molecular 
clouds or additional stars with strong winds are producing a confining ram pressure around 
this cluster, the outfiowing gas will decelerate. Measurements of this velocity will constrain 
the parameters of the surrounding material; future simulations can use these constraints to 
include the effects of this material. 



5. Discussion 

Although the diffuse X-ray flux in clusters may be used to probe stellar mass-loss, 
two notable complications in the case of the Arches and Quintuplet clusters are the X-ray 
background present in the Galactic center, and the contributions made to the overall emission 
by point (i.e., binary wind) sources. The contribution made by the X-ray background is 
difficult to quantify; Law & Yusef-Zadeh (2004) report that background contributions lead 
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to significant uncertainty in the measured X-ray flux from the Quintuplet cluster. The 
contribution from point sources is easier to handle; with its relatively high spatial resolution, 
Chandra can produce reasonable images in which the required point-source subtraction may 
be made. Point sources in the Arches cluster all exhibit a 6.7 keV line; the absence of such 
a line in the spectrum of the diffuse emission indicates that the point-source contribution to 
the diffuse emission is not likely to be significant. 

If, based on the observations made by Yuscf-Zadeh et al. (2002), we assume that Al and 
A2 are not point-like and do not subtract the contributions of point sources from the overall 
X-ray flux, it appears that the diffuse X-ray flux is consistent with higher mass-loss rates 
than the (0.3 x 10~^ Mq yr~^) 3 x 10~^ Mq yr^^ assumed for stars (less than) above 60 M©. 
As we have seen, however, the observed point-source-subtracted diffuse emission matches 
our calculated X-ray flux to within a factor 2 when we adopt the currently accepted stellar 
mass-loss rates. Indeed, lowering the mass-loss rate estimates of all wind sources by about 
30% would produce signiflcant agreement between theory and observation. But wc must 
make sure that we have a complete accounting of all the wind sources and carry through 
with a more careful point-source subtraction before we can completely conflrm such claims. 

The fact that the simulated X-ray flux density from the Arches cluster drops off more 
rapidly than the observed profile (see Figure 9) means that (1) we may have underestimated 
the contribution of the X-ray background; or, (2) we have ignored the possibly important 
dynamical influence of a conflning molecular gas outside the cluster. Interpreting data be- 
yond these radii (roughly 0.3' for Arches and 1.0' for Quintuplet) requires more detailed 
information on the cluster environment there. 

X-rays do prove to be an ideal probe of bulk mass-loss rates in clusters. The X-ray 
emission depends sensitively on the mass-loss rates of the constituent stars and, based on 
our simulations of the Arches cluster, wc can already limit the mass-loss rates to within a 
factor of 2 of the currently accepted values. This also implies that the assumed abundances in 
the Galactic center environment are essentially correct, since the mass-loss rates from stellar 
models depend on the adopted metallicity. In addition, studying the X-ray emission from 
the outer region of the cluster may eventually lead to a better understanding of the medium 
within which the cluster is embedded. This is clearly relevant to the question of how these 
unusual clusters came to be, and the relative roles played by "standard star formation" 
versus stellar capture from outside the cluster. The inferred stellar constituents of these 
clusters (Cotera et al. 1996; Figer et al. 1999a, e.g.,) also seem to be consistent with the 
required mass loss rates, so the inferred unusual mass function of the Arches, Quintuplet, 
and Central clusters continues to pose a challenge to our understanding of how these stellar 
aggregates were flrst assembled. 
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Table 1. Parameters for the Arches Cluster Wind Sources 
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Table 1 — Continued 



Star^ 




J 




V 


M 




(arcsec) 
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37 


3.54 


2.99 


5.89 


1,000 


3.0 


49 


-1.74 


14.97 


9.01 


1,000 


0.3 


61 


-1.53 


23.67 


2.02 


1,000 


0.3 


75 


7.42 


11.51 


4.86 


1,000 


0.3 


108 


7.22 


11.83 


7.08 


1,000 


0.3 


111 


0.65 


18.90 


-12.39 


1,000 


0.3 


116 


3.64 


16.48 


12.31 


1,000 


0.3 


126 


8.80 


19.13 


-0.92 


1,000 


0.3 


129 


-9.61 


10.04 


2.44 


1,000 


0.3 


132 


7.04 


20.08 


6.02 


1,000 


0.3 


149 


5.54 


21.20 


-5.47 


1,000 


0.3 


156 


5.02 


20.61 


7.93 


1,000 


0.3 



'^Numerical designations taken from Figer et al. (2002); "AR" desig- 
nations taken from Lang (2002). 

^Offset from a(2000): 17^ 45"^ 50.26^ 5(2000): -28° 49' 22".76 (Figer 
ct al. 2002). Here, positive x is ascending R.A. (to the East) and 
positive y is ascending declination (to the North). 

^Simulated with Monte Carlo. 
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Table 2. Parameters for the Quintuplet Cluster Wind Sources 
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Table 2 — Continued 



Star^ 








V 


M 




(arcsec) 


(arcsec) 


(arcsec) 


(km s~-^) 


(10-^ yr-i) 


406 


15.00 


51.50 


-26.93 


500 


0.1 



'^Numerical designations taken from Figer et aL (1999a); "QR" des- 
ignations taken from Lang (2002). 

^Offset from a(1950): 17'^43"^4.5^ 5(1950): -28° 48' 35" (based on 
Lang et al. 1999). Here, positive x is ascending R.A. (to the East) and 
positive y is ascending declination (to the North). 

'^Simulated with Monte Carlo. 
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Table 3. Properties of the Arches and Quintuplet Clusters 



Cluster 






b 




(Me) 


(Me pc-=^) 


(pc) 
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2.0 X 10^ 
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0.20 
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0.98 



^See Figer et al. (1999a). 
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Fig. 1. — Location of the 42 wind-producing stars used in the simulations of the Arches 
cluster, relative to a(2000): I?'' 45™ 50.26^ 5(2000): -28° 49' 22".76 (Figer et al. 2002). The 
radius of each circle corresponds (on a linear scale) to that star's mass loss rate. Setting the 
scale is ARl, with M = 1.7 x 10"^ Mq yr^^ 
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Fig. 2. — Location of the 31 wind-producing stars used in the simulations of the Quintuplet 
cluster, relative to a(1950): 17''43™4.5^ 5(1950): -28° 48' 35" (Lang et al. 1999). The 
radius of each circle corresponds (on a linear scale) to that star's mass loss rate. Setting the 
scale is QR2, with M = 1.5 x 10"^ M© yr^^ 
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51^ 50^ 17^ 45"^ 49^ 

RIGHT ASCENSION (J2000) 



Fig. 3. — Chandra ACIS-I 1-9 keV intensity contours overlaid on a HST NIC- 
MOS image of the Arches cluster. The image is exposure-corrected and is adaptively 
smoothed with the CIAO esmooth routine {S/N ~ 3(7). The contour levels are at 
(31, 32, 34, 38, 46, 62, 94, 158, 286, 542, 1054, 2078, 4126) x 10"=^ counts s"^ arcmin-^. 
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Fig. 4. — The 0.5-8 keV X-ray luminosity versus time from the "standard" (left) and "high- 
M" simulations of the Arches cluster. The large plot shows the variation in luminosity over 
the entire calculation, while the inset plot shows variation on a timescale of ~ 10 years. The 
winds fill the core of the cluster after ~ 2,000 years. 
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Fig. 5. — The 0.5-8 keV X-ray luminosity versus time from the "standard" (left) and "high- 
M" simulations of the Quintuplet cluster. The large plot shows the variation in luminosity 
over the entire calculation, while the inset plot shows variation on a timescale of ~ 10 years. 
The winds fill the core of the cluster after ~ 10,000 years. 
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Fig. 6. — Contours of column-integrated 0.5-8 keV X-ray luminosity per arcsec^ from the 
"standard" (left) and "high-M" simulations of the Arches cluster. In order from blue to 
cyan to green to red, the luminosities indicated by the contours are 10^°, 2.5 x 10^^, 8 x 10^°, 
and 2 x 10^^ ergs s~^ arcsec"^. The crosses indicate the positions of wind sources included 
in the simulations. 
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Fig. 7. — Contours of column-integrated 0.5-8 keV X-ray luminosity per arcsec^ from the 
"standard" (left) and "high-M" simulations of the Quintuplet cluster. In order from blue 
to cyan to green to red, the luminosities indicated by the contours are 7 x 10^^, 1.8 x 10^^, 
7 X 10^^, and 2 x 10^'' ergs s~^ arcsec"^. The crosses indicate the positions of wind sources 
included in the simulations. 
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Fig. 8. — Total 0.5-8 keV X-ray luminosity from concentric cylinders aligned with the center 
of the Arches (left) and Quintuplet clusters. The lower line in each graph represents the 
"standard" simulation of that cluster; the upper line represents the "high-M" simulation. 
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Fig. 9. — The 2-8 keV X-ray flux per arcmin^ from concentric cylinders aligned with the 
center of the Arches (left) and Quintuplet clusters. The lower line in each graph represents 
the "standard" simulation of that cluster; the upper line represents the "high-M" simula- 
tion. Crosses and error bars in the left graph represent flux measurements from Chandra 
observations of the Arches cluster. 
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Fig. 10. — Average radial velocity of gas in the "standard" simulations of the Arches (left) 
and Quintuplet clusters, versus distance from the center of each cluster. 



